Shape-invariant quantum Hamiltonian with position-dependent 
effective mass through second order supersymmetry 



A Gangulyi fl and L M NietoH 

^Departamento de Fi'sica Teorica, Atomica y Optica, Univ. de Valladolid, 47071 Valladolid, Spain 
2city College (C. C. C. B. A.), University of Calcutta, 13 Surya Sen Street, Kolkata-700012. India 

Abstract. Second order supersymmetric approach is taken to the system describing motion of a 
quantum particle in a potential endowed with position-dependent effective mass. It is shown that 
the intertwining relations between second order partner Hamiltonians may be exploited to obtain a 
simple shape-invariant condition. Indeed a novel relation between potential and mass functions is 
derived, which leads to a class of exactly solvable model. As an illustration of our procedure, two 
examples are given for which one obtains whole spectra algebraically. Both shape-invariant potentials 
exhibit harmonic-oscillator-like or singular-oscillator- like spectra depending on the values of the shape- 
invariant parameter. 



PACS numbers: 03.65.Ca, 03.65.Ge 

1. Introduction 

The effective mass (EM) Schrodinger equation with position dependent mass is a very useful model 
in many applied branches of modern physics, e. g. semiconductors (T], quantum dots |2| , and '^He 
clusters [3J. In these cases, the envelope wave function actually provides a macroscopic description 
of the motion of carrier electrons with position dependent (or equivalently material-composition 
dependent) mass. Recent interest in this field [4-17] stems from extraordinary development of 
nanostructure technology. This sophisticated technology of semiconductor growth like molecular 
beam epitaxy technique, makes production of ultrathin (nonuniform) semiconductor specimen a 
reality nowadays. Consequently the study of such EM equation becomes relevant for deeper 
understanding on the non-trivial quantum effects observed in those nanostructures. Another area, 
where this equation has extensive use, is the study of quantum many-body systems in solids. For 
instance, the computational limitation of the Green's function quantum Monte Carlo method is 
removed by using the so-called pseudopotentials. The EM equation appears in the process when one 
attempts to replace the nonlocal (and thus disadvantageous) angular- momentum projection operator 
of the pseudopotentials [TSl . 

On the other hand, supersymmetric (SUSY) approach ||T9l has been proved as a powerful tool 
in quantum mechanics identifying the isospectral partners in bosonic and fermionic sectors, and 
thereby generating a hiererchy of solvable Hamiltonians f20l. In the standard SUSY approach for 
EM Hamiltonians the ladder operators are taken as first order derivative operators similar to constant 
mass (CM) case, but now they depend on both superpotential and mass function. As a result one 
obtains two partner EM potentials with the same effective mass sharing identical spectra up to the 
zero mode of supercharge. In the literature a generalization of standard SUSY known as higher 
derivative SUSY (HSUS Y) or ?^-fold SUSY, find a good usage for CM models [21-31]. This method 
keeps the basic superalgebra and differs from standard first order SUSY in that the supercharges are 
represented as "N-th order (N > 1) differential operators. Recently this higher order approach has 
been extended for EM Hamiltonian llT6l . When one is interested in the solvability of the Hamiltonian, 
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shape-invariance (SI) is an important criteria |[32| , because in this case one can obtain full spectra by 
successive application of raising operator. For instance, some classes of SI Hamiltonian have been 
discovered recently through first order SUSY formalism 161 [T2l [131 . 

Our purpose in the present work is to extend the concept of SI to second order SUSY (SSUSY) 
in the effective mass framework. This type of extension has been done very recently ll26l for CM 
case. We have derived a new relation between the potential and mass functions, which gives an 
exactly-solvable (ES) model satisfying a simple SI condition. Two specific examples are given, one 
is hyperbolic and other is algebraic. Both these Hamiltonians are new ES models, because their wave 
functions and energy eigenvalues can be obtained by purely algebraic means. For certain values 
of the characterizing parameter these potentials acquire an inverse-square singularity at the center. 
Such type of singular potentials has important applications in various fields e. g. molecular and high 
energy nuclear physics 1 33 1 . In the constant mass limit both of these reduce to well-known harmonic 
oscillator or singular oscillator depending on the values of the SI parameter. 

In Section |2] we set up first order SUSY for EM models. Section [3] is devoted to SSUSY. 
After getting compact expressions for intertwined Hamiltonians in terms of superpotentials and mass 
function, we obtain the solution of zero mode equation of supercharges. In Section |4] we study in 
detail the connection of the second order scheme with the standard one. Section |5] contains the 
main result concerning SI for SSUSY scheme. A new relation between potential and mass function 
leading to SI model is derived. Two examples satisfying such relation are provided. Finally, Section|6] 
contains the concluding remarks. 



2. First order SUSY and EM models 



The well-known superalgebra defined by the super Hamiltonian Hg and the supercharges Q, js 

H,^{Q,Q^}, Q-" = Q^" = [Q,Hs\ = [Q^,H,]^0, (I) 

where the symbol 'f ' denotes the usual hermitian conjugation. The supercharges are represented in 
one dimensional quantum mechanics by the following 2x2 matrices 

At 




A 







Consequently, the super Hamiltonian Hg is diagonalized as 



Hg 









i7_ 



PJK i?- = AA 



(2) 



(3) 



In the standard SUSY approach to EM models the ladder operators A, A^ are first order differential 
operators [jjl 



1 



5 + W(x), 



1 



d + W{x) 



1 



\/m{x) 



dx 



(4) 



m{x) ^ ' yjm{x) 

where W(x) is the EM superpotential and the prime denotes derivative with respect to x. Note that 
we have omitted the factor fi^ by defining the atomic units such that — 2. Denoting H± as 

1 



-d 



m{x) 



d]+V±{x), 



(5) 



the realization Q straightforwardly expresses the Schrodinger EM potentials in terms of mass and 
EM superpotential 



V+{x) = - 



W(a;) 
VMx) 



V-^{x) = V+ix)+2 



W'{x) 1 
y/mix) y/m{x) 



1 



(6) 



The representation of ladder operators is not unique for EM models, for other variants see I12II13I . 
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The main feature of above construction is the intertwining relations between the bosonic and 
fermionic parts, which plays a crucial role in SUSY theory 

AiJ+ = i/_A, iJ+At = AtiJ_ . (7) 

The vacuum states may be characterized as 

AV'+(a;) = 0, or A^^ (a;) = 0, (8) 

where 'ip^ , n 0, 1 , . . . denote the bound state wave functions for the Schrodinger Hamiltonian H± . 
Let us point out that physically the mass function m{x) is positive definite and finite everywhere in 
the domain of definition of the EM Schrodinger equation. Now, the existence of a zero mode will 
depend on the asymptotic nature of the superpotential W(a;) and mass function m{x). The ground 
state wave functions may be computed from equation dHJ 



iPq{x) oc exp 



y/m{T)W{T)dT 



ipQ {x) oc -\/ m{x) exp 



^m(r)W(r)dr 



(9) 



The existence of the vacuum states is ensured if and only if the integrals satisfy the following 
asymptotic condition 



\/ m{T)W{T)dT > + oo ( or — oo) . 

X — >±oo 

Clearly both of the above conditions can not exist simultaneously and so at most one of the zero 
modes ([9]) will be annihilated in which case SUSY is unbroken (or exact). On the other hand, if 
neither of the zero modes exist then SUSY is called spontaneously broken. In this context it should 
be kept in mind that this argument fails for periodic Hamiltonians ll34l because in that case one 
usually considers Bloch solutions and so the square-integrability criteria disappears [35-37]. In fact 
the question of existence of zero modes are related with the specified functional space on which 
a quantum system is to be considered |30| . However in this treatise we will not consider periodic 
models. Thus for the systems defined on linear space it must be safely concluded that two zero modes 
can not exist simultaneously in the first order SUSY formalism. In the next section we will show that 
this is not the case for the SSUSY scheme. 

For definiteness, let us suppose that the bosonic sector (i?+) is fully known and possesses 
normalizable zero-energy state. Thus, if either ^p'^ or W{x) is known, the other can be obtained 
exactly by formulae (|9|l. Now, one can extract full knowledge about the fermionic sector (i?-) by 
the use of the intertwining relations (|7]i: 

(x) = A^++i (x) , /b^F {x) = At (x) ■ E- = E+^, ,E+=0, (10) 

for n = 0,1,2.... 

Let us now consider the EM eigenvalue equation 

H^u{x)i}{x) EE [T^u{x) + V^m{x)]'P(x) = eV'(x) (11) 

The first step to study this equation is certainly to choose a suitable form of the hermitian kinetic 
energy operator T-^m {x). There is an intrinsic ambiguity in selecting such form as this class of physical 
problems are suffered from non-commutativity of momentum operator p = —i\/2d and the effective 
mass operator m{x). Several forms had been proposed in the literature for T^m{x), and considerable 
efforts were made to remove the non-uniqueness of the kinetic energy operator [38-43]. But still the 
problem of ambiguity remains a open question in this field. However one may rely on almost general 
representation suggested in Ref. |39| 

Tem(x) = — (jn'^p m}^ p + m^pm^pTn^^ , a + b + c = —1 . (12) 

The parameters a,b,c in the above equation are usually called 'ordering parameters'. Most of the 
present authors prefer to start from (fT2] | as this representation includes many special forms used in 
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different context. For instance, the authors in Ref. ||6l considered following kinetic energy operator 
to apply first order SUSY 

2 \m J 

which was first proposed by BenDaniel and Duke ||38l . and is contained in the general representation 
(fT2] ) for the special choice a = c = 0,b = — 1. In this connection we would like to mention an 
interesting work |12| where the authors have proposed a new representation of first order SUSY 
ladder operators including ambiguity parameters of the kinetic energy operator ( fT2l i and have 
obtained a substantial generahzation over the result of Ref. 161. It is now well-known that lfT3HT4l the 
representation of ladder operators can be made free from the ambiguity parameters [see equation ©I 
at the cost of constraining the local potential Vem{x) in equation (fTTT i with a pseudo potential term 
p{m) thereby considering the so-called effective potential energy VEa{x). For the two-parametric 
representation (fT2] | of the kinetic energy operator this pseudo potential term is given by 

1 + b m" m'^ 

P{m) = ^ 2"^^' ?? = l+b + a(a + b + l), (13) 

where we have absorbed the parameter c using the constraint a + b + c = —1. In what follows we 
will consider the following general EM Schrodinger equation 



m{x) 



V'(x) = eip{x) , V^^ix) = V^„{x) + p{m) . (14) 



In above equation, Km (a;) represents the local potential strength and the pseudo potential p{m) is 
given by (fT3] l, the latter depends on the ordering parameters a, b. We will assume that a,b are real. 
The EM Hamiltonian (fT4l l may then be identified with bosonic partner Hamiltonian iJ+ in (|5]l as 

V+{x) = V^»{x) - e , e<E+=0, (15) 

by assigning Vem{x) = V+{x) — p{m) + e. The quantity e is often termed as factorization energy in 
the SUSY procedure 1441 . 



3. Second order SUSY and zero mode equation 

We will now replace the intertwining operators A, in equations ([T])-© by the following second 
order operators 



- +W{x)d + c(x) , ^-d^ 
m m 



W (x) + 2 



9 

m 



d+[c{x)-W'^{x)] , (16) 



Wm{x)=V^{x) + ^-^. (17) 



in which we have used the abbreviation 

m'{x) 
rTn?{x) 

In the following we will use the terminology "superpotential" for the function Wm(a;). Clearly the 
super Hamiltonian Hs given by ([T])-©, is now a fourth order differential operator, and it will be 
physically meaningful if it can be expressed as a quadratic polynomial in the usual EM Schrodinger 
Hamiltonians. Thus we will introduce 

"o" L)^ = -a (la) (18) 
Our task is to determine the following matrix identity 

Hs = hi + hh, + hh , (19) 

where li, I2 are arbitrary fixed real numbers. After some involved but straightforward steps one may 
express v±{x) in terms of c{x), the superpotential rn{x) and the mass function m^x) 

v±{x) = {\t T + - c{x) ~ I , (20) 

2 2m 2 2 
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where the function c{x) is given by 

^ ~ 2 4 2mW™ TO V2W„J 47n3 2to2 m V2W, ' • ^ 

In above equations the quantity K^^ is taken, without loss of generality as 



K:, = +^ll-Ah. (22) 

It should be emphasized that the quantity may be purely real or purely imaginary according as 
l\ > 4/2 or ll < 4/2- It is that constant which plays the role of determining reducible or irreducible 
SSUSY II22I . In the next section we will show that forhermitian quantum mechanics, SSUSY scheme 
could be reduced to first order SUSY for real K3 only. 

Note that the EM Hamiltonian hg always commutes with super Hamiltonian Hg due to 
the relation ( fT9l l. and hence both have simultaneous eigenstates. Thus the intertwining relation 
(|7]l implies that Schrodinger Hamiltonian (fTSl ) is doubly degenerated (up to zero modes of the 
supercharges), and wave functions are connected according to (fTOl i. where now denote 
eigenvalues of h±. It should be mentioned that above expressions for v± and c{x) satisfy additional 
intertwining relations between EM Hamiltonians h± 

Ah+ = h-A , h+A'^ = A^'h- . (23) 

This relation is important not only because of its elegant description, but also one can start from the 
requirement ( |23] |. and may obtain expressions ( |20| | and (l2Tl l. 

We will now show that zero modes of both operators ( fTSI l may exist simultaneously, in contrast 
to the standard first-order SUSY. To understand this let us write down the normalizable solutions 
of zero mode equations (O for supercharges. Note that the equation ^ becomes a second order 
differential equation with the replacement of A, by A, A'^ of (fT6] i. These equations may be brought 
to the form similar to CM Schrodinger equation 



2W„ J ^\ 2W™ 
by the transformations 

V'^(x) = y^Mxj^ix) elTir ™MW,„(r)rfr] (25) 

The normalizability of tp^ clearly depends on functional forms of both superpotential W ,„ (x) and 
mass function m{x). To illustrate it, consider that (t){x) in (|24] | is identified with normalizable ground 
state wave function of a known solvable CM Hamiltonian by suitably choosing the superpotential 
mix). Then the prefactor (t>{x) of V'o^ in ( l25T l is well-behaved, and hence both if)^ will be 
normalizable if m{x) is also well-behaved and the integral J toW mdr is finite. Thus in this instance 
both operators ( fT6b have normalizable zero modes, as was the situation for CM models. The explicit 
solutions for zero-mode states of both supercharges may be written as follows 

V.±,.(x)ocy;^e[^r^^M'^^ J = l,2. (26) 

2 W rn 

Note that for = 0, Fi (x) = F2 (x) and so in that case both operators may have at most 
one zero mode. It should be mentioned that both of the zero modes are also formal eigenstates of 
Schrodinger Hamiltonians h± for real A'3 only : 

h±i^o,j = VjAj > Vj = ^ • (27) 

Before concluding the section it may be pointed out that the quantum systems built upon 
K-th order representation of ladder operators were categorized as type A K-fold SUSY (see for 
details Refs. Il30l [T6l ). Hence it is not difficult to show that the systems investigated in this 
article are generically special cases of type A 2-fold SUSY. Furthermore in Ref. (30\ it was shown 
systematically that the zero modes of one higher order supercharge can admit both physical and non- 
physical states, the latter may be used as a good transformation function to develop new solvable 
system 1371 . 



(j){x) = , (24) 
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In the previous section we have derived the expressions for second-order SUSY partner potentials 
v±. In this section we want to investigate whether v± may be expressed through standard SUSY 
formalism. That is, to say whether there exist superpotentials Wi, W2 in terms of which v± may be 
brought to the first order form given by (|6]l. To proceed systematically, let us write h± as 

Ki, h+ = h^_^^ = A^Aa + K2, Aj = + , (28) 



,(1) 



AiA 



where A^ are first order EM ladder operators and Ki , K2 are suitable constants to be determined. In 

('2) 

h\_ interchangeably for 2-SUS Y partners 



this section we will use the notations 



and h+ 



h±. Let us consider that K3 is purely real. Assuming for definiteness K3 > 0, two types of solutions 
for Wi, W2 are possible. 



(a) Type I: Ki 



K3+I1 



K3-I1 



2 ' 2 
The first order superpotentials are given by 



Wi, = 



± 



2W, 



(29) 



The relation between second order operators A, and the first order operators Ai, A2 in this 
case is 

yiU = (AtA2)(A^A2 
(b) Type II: K2 = Ki = 

K3 



K3), AA^ = (AiAI)(AiA1 - Ks). 



(30) 



Wi,2 = 



1 



2W, 



± 



2W, 



TO 



1 



(31) 



The operator relation (l30l ) becomes 

yiU = (AtA2)(AtA2-i^3) 



AA'' = (AiA1)(AiA1 -iiTa 



We see that both types of reductions are distinct for the superpotential W2 unless K3 = 0. One 
may check readily that Type I reduction allows factorization of ladder operator A in terms of first 
order operators Ai , A2 as 

(follows from the identity: a| Ai 



A A1A2 



A2AT 



However this is not possible for Type II reduction. But in both cases v± are in 1-SUSY form, which 
is actually important. Their respective 1-SUSY partners may be written down using the formula (|6|l 



„(i) 



Wi 

TO 



Kl: 



(2) 



2_ 

TO 



1 



1 



TO 
(1) 



1 

,(2) 



+ K2.{32) 



According to standard SUSY, and are isospectral where v_ = w_ and v_^_' = f + are given 
by (|20] |. Hence we see that second order SUSY formalism may give us opportunity of studying two 
standard SUSY pairs (/I'f''', /I'f '') simultaneously for j — 1,2. The ground states of both pairs of 
Hamiltonian are zero modes of first order operators Aj , Aj given by ( |28] l. These can be computed by 
substituting for Wi, W2 from ( |29] ) and dSTT i (corresponding to Type I and Type II reduction) into the 
formulae (|9|l for W. For instance the zero modes for the opertors Aj,Aj,j = 1,2, corresponding to 
Type II reduction are 

-r F2(r)dr 



'(1) 



'(2) 



OC 



(33) 



and 



where V'o 2 ^i"^ given by 

(i) (i) 

Hamiltonians ,h_ for j = 1, 2 respectively. 



-U) 



denote ground state wave functions for the 
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Let us point out that for imaginary K^, given by (l22T i. the above reduction is not possible in 
hermitian quantum mechanics ||3T1 . because both first order superpotentials Wi, W2 will be complex 
[see equations (|29] l and (|3T|i1. In fact this will lead us to an irreducible transformation between real 
and complex potentials for EM Hamltonians ||45| . 

The discussion in this section clearly shows that the factorization of a higher order linear 
differential operator in terms of lower order operators is in general non-unique. That is to say, one 
operator can admit both reducible and irreducible representations. Hence from strictly mathematical 
sense the concept of reducibility of HSUSY must be defined on the basis of an additional restriction 
on Hamiltonians that they are factorizable according to ( |29] l or dSTI ). In the next section we will 
propose a higher order SI criteria for EM Hamiltonians. Our purpose of introducing this section is 
to compare the result of SI obtained through HSUSY with that obtained via first order SUSY. This 
will give us a better insight about why and how SI scheme proposed in this article is an important 
generalization over the SI formalism in the standard approach. 



5. Higher order SI criteria for EM Hamiltonian 

In Ref. |l6l two types of SI criteria were discussed for EM Hamiltonian generated through first order 
SUSY let apart the generalized treatment proposed in Ref. {lT\. Recently a kind of deformed SI 
criteria has been introduced |,13J, which is also in the context of standard SUSY. Searching for SI 
Hamiltonian is useful, because this integrability condition leads to certain relation between potential 
and mass functions producing an ES model. For instance, claiming first order SUSY partners 
y_(a;; A) — Vj^{x] A) = 2A, from equation (|6]l, one obtains following relation between first order 
superpotential W(a;) and mass m{x) : 

W(a;) = i (^-^^ +\j V^dT. 

Just this relation was discovered in Ref |6|. Here we wish to enquire this simple SI condition for 
SSUSY pai-tners w-), given by (l20l i. 



5.7. Theoretical construction 
Let us consider 

h^{x;X) ^ h+{x;X) + 2X, A>0. (34) 

It is not very difficult to see that the above requirement expresses second order superpotential in the 
form 

W™(:r) = 4=> 9{x)=j + X r V^dT, (35) 
y'm(x) J 

where 7 is an integration constant. Substituting this expression for W into ( |20] | for v+{x), we 
obtain following relation between potential and mass 

v+{x;X)^^ r^ + ^-T7^-U + l7 ' (36) 

where we are considering K3 > 0. 

We stress that this is a new relation between potential and mass which gives us a class of SI 
Hamiltonian h+. The whole set {?/'+(a;; A), E"^} of eigenstates and spectra for /i+ can be constructed 
by exploiting the intertwining relation (|23] l. The procedure is, as in the case for harmonic oscillator, 
to apply successively the raising operator A'' upon zero mode of the lowering operator A. Note that 
the second order operator A has two zero modes ip'^^ , j — 1, 2, given by ( |26] l. Both these zero 
modes are formal solutions of Schrodinger equation for h+ [see equation ( |27] |1 with the eigenvalues 
771 = (— /i + K^)/2 and 772 = — (/i + A'3)/2 , where 771 > 772. Hence we obtain double sequences 



,71 = 0,1,... (38) 
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of eigenstates based on both zero modes. The labels of eigenstates will depend on the values of the 
SI parameter A in equation (|34] i. For simplicity let us first consider the case 

A > , i^s > (37) 

In this case the wave functions and energy eigenvalues of h+{x; A) may be expressed as follows 

i>tni^;X) = (yit)"^+2(a;), ^m + 2nX., 

i^tn+i{x-A) = (yit)"VM(a:), ^2^1 =77i + 2nA, 

We will now turn to the case for A < K3/2. Note that for all values of A(> 0), the ground state 
will be given by the zero mode ipQ2 with the eigenvalue 772- Higher excited states will be obtained 
in a similar procedure as described above, but one has to relabel the states according to the range of 
values of A. To illustrate let us take the most general situation where 

K3/2{k + 1) < a < K3/2K , /t = 1, 2, . . . (39) 

Here first k members will be given by the sequence based on the zero mode ipQ2- Thus the wave 
functions are 

V'+(a;;A)-| (yit)'^7/-+2(a;), n = K + 2,K + 4,... (40) 

with the spectra 

{rj2 + 2nX , 7T, = 0, 1, . . . K 

7?2 + (n + K)A, 71 = K + 2,K + 4, ... (41) 

rji + {n — K — 1)A , n = K + 1, K + 3, . . . 

It is a trivial exercise to verify that for A = K3/2K, the sequence built on ipQ^ coincides with that 
built on ipQ2 except for the first k members, which are the lowest states and only singular solutions. 
The wave functions and spectra will be provided by single sequence 

A=^: i>t{x;X)^{A^)"ij+2{x), E+ ^ ^ + , n ^ 0,1, . . . . (42) 

It should also be kept in mind that for SSUSY scheme, zero modes of both operators A and 
may simultaneously exist. This implies that the above sequences may terminate at n = N{> 1) if 
{A'')^'ipQj oc tp^ji , j, j' — 1, 2, tpQji being the possible zero modes of A^, given by ( l26b . Note that 
for A = K3, we have from (l38b harmonic-oscillator-like spectra 



^n = -(A + M + U+i) A, n = 0,l,. 



while for A 7^ K3, the spectra given by ( l38T l, i4T[ or ( l42b separately coincide with those of singular 
oscillator. 

Under the SI condition ([34]), all levels can be obtained in a closed analytic form. For instance, 
first three members of the sequence { (A^) 'ipQ2 i^)} 



y^+2{x) « m^/Hx) [g{x)]^ e^^ , yLt7/;+2(x) cx [g^ix) + K3 - 2X] i,+2i^) , 

(yit)' Vo+2(^) « ligH^) + K3-Axy + 2 a {k^ - ax)] i^+2{x) . 



(43) 



Higher members can be constructed in a similar fashion. The members of the other sequence 
{ (A^)" rpQi {x)} can be obtained by simply changing the sign of K^, in the corresponding members 
of the sequence {(A^^" ■iI}q2{^)}- The nature of the function g{x) given by (l35T l is crucial for 
the normalizability of the wave functions. In the first place, ■m{x) must be so chosen that it 
verifies g'^{x) +00 as a; — > ±00 and remains finite otherwise. Secondly, suppose that g{x) 
is nodeless in the whole domain. Then the potential v+{x; A) is non-singular and both sequences 
provide regular (non-singular) solutions. But if g{x) has a node at x — xq then the potential is 
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singular at that point for A 7^ K^. Let us note that the node of g{x) is of first order, because 
g'{x) = m{x) is always nodeless for chosen mass function. Hence near x ^ x^, v+{x) will 
behave like C(x— a;o)~^, where the constant C characterizes the strength of the singularity. It is well- 
known that self-adjoint extension of such Hamiltonian can be determined on the whole domain for 
the range — 1/4 < C < 3/4. The singularity will be attractive or repulsive according as A > K3 > 
or < A < ifa. One may readily check from (l43T l that V'o^2(^) singular for A < K^, and so in this 
case depending on the values of A and other parameters some or all of the members of the sequence 
{ ijjo2 i^)} have to be deleted from the set of regular solutions of h+ {x; A). 

Hence we have proved that the Hamiltonian A) with the potential v+{x; A) in (l36T l is a 

new ES model possessing second order SI condition in the EM framework. Two remarks are in 
order. As one lowers the value of A by increasing the value of k in (3% . the levels become closer 



and closer and in the limit A 0{k — + 00) only two levels will be left namely 'ipQ2{^) ^pQ^{x) 
with the eigenvalues 772 , rji respectively. Clearly this gives a quasi-exactly solvable system with two 
known levels. Since in this article we are only interested about ES models, k will be a finite quantity 
or equivalently the SI parameter A is a strictly non-zero finite positive number Secondly we have 
already mentioned that if the function g{x) in ( |35] | has a node at a; = xq, the potential u+(a;; A) given 
by (|36] | will be singular at that point. Hence extreme care should be taken to decrease (or increase) 
the value of A in order to keep the strength C of the singularity controlled (—1/4 < C < 3/4), as 
for a given mass function m{x) and the constant the strength C is inversely proportional with A. 
This will be clear if one expands the function g{x) about its node xq giving 

C=(^^/^. (44) 



4m(.To) 

Next we are going to construct two classes of examples based on the theoretical model proposed in 
this subsection, where this last comment will be crucial to get the physically acceptable Hamiltonians. 



5.2. Hyperbolic mass and potential 
Let us choose the mass function as 

m(x) = (a + /3tanha;)^ /3 7^ 0, |a| > |/3| . (45) 

This mass depicts a smooth step function. It increases (or decreases) from the value [a — /3)^ for 
X = —00 to the value {a + /3)^ for x — +00 according as a/3 > (or aP < 0). From dSST l the 
function g{x) may be computed 

g{x;\) = 7 + A [ax + /? In cosh x] . (46) 

The potential turns out to be 

-2 ~ li\ ^ /3sech^a:;[/3tanh^ X - 4atanhx - 5/?] 



v+ix;X) = — 



^+2 



4(q; -I- /3tanhx)^ 



(47) 




V, (X) 
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Figure 1. The effective mass J45b and potential function ( 147 1 with Zi=5,/3 = 7 = 1,A = 
= 4 that depicts harmonic oscillator-like barrier. The graphs are drawn in the atomic 
units (a.u.) defined by fi,^ = 2 for four different values of the parameter a. 
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The nonsingular SI potential ( |47l ) and effective mass ( |45] ) are shown in Fig. [T| corresponding to 
Ql) [3 — ^ — 1, Zi = 5, A = Kj, ~ 4 and (ii) four different values of the parameter a. For all 
values of a (with the restriction \a\ > \(3\, (3 0, which is necessary to have nodeless nonconstant 
mass) we have a harmonic oscillator-like well. As the magnitude of a decreases the well becomes 
flatter. The constant mass system will be recovered with (3 = 0, and consequently the potential 
will reduce to harmonic oscillator for A = K3. For very small values of f3 the wave functions 
tend to those associated with harmonic oscillator. On the contrary, for large values of (3 the wave 
functions tend to spread over the region whose size grows as (3 increases due to the influence of 
step mass (l45T l in the potential profile. This behaviour is demonstrated in Fig.|2] which contains the 
potential and wave functions associated with ground and first two excited states corresponding to (i) 
a = 7 = 1, A = isTa = 4 and (ii) /3 = 0.9, 0.001. 



i/fn (X) 




Figure 2. Wave functions associated with ground and first two excited states and the SI 
potential ( 147 1 corresponding to /3 = 0.9 and /3 = 0.001 while other parameters are 
Q = 7 = 1, A = K3 = 4 (a.u. are used). 

The nature of the potential is dramatically changed for A 7^ K^. As can be seen from ( |46] l that 
near x ~ xq, the potential behaves like 

A) ~ - — ^—^ , 

[x - XqY 

where x = xq \s node of g{x) and the strength C of the singularity is to be computed from the 
formula ( l44b . For the weakly attractive singularity (—1/4 < C < 0) the physical wave functions 
tend to vanish atx ^ xq according to the boundary conditions. This behaviour is clearly appreciated 
from Fig. |3] where we have plotted singular and non-singular potentials with their probability density 
functions (|-0„(a;)p) associated with first three lowest levels corresponding to (i)A = 4and (ii) A = 6. 
The other constants are taken as li = 5, K3 = 4, /? = 7 = 1, a = 2, which gives C = —0.04 and 
xo = -0.09. 




Figure 3. Non-singular and singular (attractive) potential 147 1 along with probability-density functions 
(in a.u.) associated with ground and first two excited states corresponding to A = 4 and A = 6 for 
= 5, /sTa = 4, a = 2, /3 = 7 = 1. 
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5.3. Algebraic mass and potential 

As our second example, let us take the following mass function studied in Ref. 161 

m{x)^[^^^ ' ">0' (48) 

This mass function remains strictly positive definite everywhere and approach a constant value 1 at 
both infinity. The SI potential is 

where 

A) = 7 + A[a:: + (a - l)tan"^a;] . (50) 




Figure 4. Non-singular potential ( 1491 ) and mass function ( 148b . depicted in a.u., for l\ = 
5, 7 = 1, A = i^s = 4 corresponding to (i) q = 0.3, 0.6 (first row) and (ii) a = 1, 10 
(second row). 

The non-singular potential is given by ( |49] l for A = K-^. For < a < 1, it is a bistable 
potential and for a > 1 it is a single potential well, while a = 1 recovers the constant-mass system 
giving harmonic oscillator well. Fig.|4]describes non-singular potential (|49] l and mass function (|48] | 
corresponding to (i) a = 0.3, 0.6 (first row) and (ii) a = 1, 10 (second row). From this Fig. it is 
clear that for small values of a, the height of separator between two wells grows producing a thin 
barrier. On the other hand for very large values of a, the single well becomes sharper compare to 
that of standard harmonic oscillator well. 

The potential acquires an inverse-square singularity for A 7^ K^, which is repulsive (attractive) 
for < A < K-i (A > K^). For weakly repulsive singularity (0 < C < 3/4) the regular (non- 
singular) wave functions are given by the sequence {(^^) ^^0^1(2^)}- According to the boundary 
conditions required to have self-adjoint extensions of the Hamiltonian, these wave functions tend to 
vanish at x — > xq, xq being node of the function g{x) given by ( |50] l. This behaviour is depicted in 
Fig.|5] wherein we have plotted singular potential along with probability density functions associated 
with first three excited levels of the sequence {ipn{x; A)} given by (|42] | for A = 2, other parameters 
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being taken as /i = 5, Ky, ~ A, a = 2,7 = 1. These choice of parameters yields the strength C and 
center Xq of the singularity as C = 0.25, = —0.26. 



v,(x) |(i/n(x) 1-= (n = l,2,3) 



-2 ^ -l 


'P 




.0.6 


--l^il^ 

|l/'2l' 




. 2 




-0.26 2 



X 



Figure 5. Singular (repulsive) potential l |49t and the probability-density functions (in a.u.) 
corresponding to non-singular wave functions associated with first three excited levels for 

/i = 5,a = 2,7 = 1,A = 2,JC3 = 4. 

Physically the well confining the particle is divided into two zones by a high thin barrier at 
Xq = —0.26. But note that for weakly repulsive case this barrier is penetrable by a quantum- 
mechanical particle. 



5.4. SI condition in reduced first order SUSY 

In subsection lS.ll we have proposed a second order SI condition ( l34b and consequently have obtained 
a new system of ES Hamiltonians h+{x; A) where the potential v+{x; A) and the mass m{x) are 
related according to ( |36] |. Here we wish to show that the same system ( |36] | (up to a constant 
shift) could be studied through first order SI formalism provided the Hamiltonians are factorizable 
according to ( |29l ) or dJTT ). To see this, consider the 1-SUSY pair (i^-i- = ^ , ^ ) given by < l201 t and 
(3% corresponding to Type II reduction ( |3TI ). Here the SI parameter will be K3 instead of A, while 
the other parameters have to be considered fixed. For convenience, in the following we will use the 
symbol ao for K3. The relation between first order superpotential W2(x; ao) and the mass function 
m{x) reads from dsTl ) and dSSl l 

To have zero energy ground state, consider the Hamiltonian lv_^' {x; ao) = h^_^_ {x; ao) — ^^2(00) = 
a\{x; ao)A2{x; ao), where the operator A2 is given by ( |28]) and K2{ao) = — (ao -|- ?i)/2 for Type II 
reduction. Let us use the following abbreviations for the shifted potentials 

v':^\x;ao) = i;*^^-'(a;;ao) - K2iaQ) , uL^^(x;ao) = v^_^\x-ao) - K2{ao) , 

where (a;; ao), en correspond wave functions and eigenvalues. Note that the wave functions 
and energy eigenvalues of h^{x] A), given by ( fTSl l and ( [36l ) are related with ao), e,t'"^''} by 

^+(a;;A)cx0+(2)(x;ao), i?+ = e+(2) + , n = 0, 1, . . . 
where ao = K^^. One then find from (ISTT i following SI condition in the reduced first order SUSY 

(2) f2l 

u_ "^(a;; ao) = (x; ai) + _R(ai) , ai = /(ao) = 2A — oq , -R(ai) = 2A — ai . (52) 

~ (2) 

Thus full spectra of h\_ {x; ao) can be recovered according to standard prescription 



"=0 



-(2) _ n .+(2) 



0, £« = I]fe=i^("fc)' ^fc /(flfc-i) 



(?!);l;^^^(a:;ao) cx A^(a;; ao)Aj(x; ai) • • • aJ(x; a„_i)0o (a;; a„) 



1,2. 
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where (a;; oq) is zero mode of the operator A2(x; ag), given by (|33] l. Two points are to be noted. 
Firstly imposing the restriction of factorizabihty on the Hamiltonians one would come out with the 
SI condition ( |52] | and the relation ( BTT l. which have not been reported so far in the literature for EM 
Hamiltonian. Secondly the proposed SSUSY scheme definitely leads to a generalized SI criteria ( |34] |. 
because in this case no such restrictions need to be irnposed on the Hamiltonians. 

Our final remark is that in the CM limit m llj], both potentials A), given by ( |47] | and 

(|49] ) reduce to known SI potential v™{x\ A) : 

.+ (.,A) = ^ ____(^A + -j , (53) 

and the wave functions and energy eigenvalues exactly coincide with those B6ll47]| for i53[ . 



6. Conclusion 



In this article we have used SSUSY scheme for describing dynamics of a quantum particle with a 
position-dependent mass. We have derived a compact expression of 2-SUS Y pairs {v+ ,v^) in terms 
of second order superpotential Wm(x) and the mass function m{x). A detailed analysis has been 
given about zero mode equations of second order supercharges and possible reduction of SSUSY 
scheme to first order SUSY. Based on the existence of intertwining relation between 2-SUSY partner 
Hamiltonians h±, we have obtained a new relation between potential v+{x) and mass m{x) leading 
to a simple SI condition. As a result full spectra is achieved by successive application of second order 
raising operator A"^ upon the zero modes of the lowering operator A. It is shown that in the reduced 
first order SUSY approach, one may obtain a new relation between first order superpotential W2 (x) 
and mass m{x) and the same system (up to a constant shift) could possess a first order SI property 
provided the Hamiltonians are factorizable. The advantage of using SSUSY scheme is that it deals 
with a generalized but simpler SI requirement, namely partner Hamiltonians differ by a constant A. 

We have constructed explicit examples for two types of position-dependence of mass function, 
one is hyperbolic and the other is algebraic. The corresponding potentials have very different aspects 
based on the values of the SI parameter A. Both are non-singular for A — K^, where the constant K3 
appears in the process of relating Schrodinger Hamiltonian with super Hamiltonian. The parameters 
f3 and a characterizes position-dependence of mass ( |45] | and (08]) respectively. The non-singular 
algebraic potential ( |49] | shares same qualitative features as those discovered in Ref. ||6l. For A 7^ K3, 
the potentials acquire an inverse square singularity, which is attractive (repulsive) for A > K3 
(A < The week strength of singularity (—1/4 < C < 3/4) is particularly of physical interest. 
Because for weakly attractive singularity (—1/4 < C < 0) the ground state energy remains finite. 
On the other hand, for weakly repulsive singularity (0 < C < 3/4) the barrier is penetrable by 
quantum-mechanical particle. In both instances, regular (non-singular) wave functions vanish at the 
center of the singularity. 

It is important to clarify that the novelty of our work lies in the fact that we have generalized 
for the first time the concept of SI to the second order SUSY approach for EM Hamiltonians. In 
this context one should look into the explicit classification done in Ref. Ill 61 via INf-fold SUSY 
approach to EM quantum systems. The systems investigated in this article belong to a special 
subclass of the systems studied there, which possess second order shape-invariance. Hence one has 
a significant advantage of obtaining wave functions and energy eigenvalues for such Hamiltonians 
over the existing general method for arbitrary potential and mass functions. 

We would like to mention that in many practical applications, where continuous spectra is of 
interest E], SSUSY scheme used here may be utilized to relate transmission and reflection amplitudes 
of partner Hamiltonians. The idea of SI condition for SSUSY scheme can be straightforwardly 
extended to general 3\f-th order representation (N > 2) of ladder operators. However obtaining a 
physical model will be much more difficult, because one has to make an appropriate choice of the 

K m ^ 1 is taken to tally with the conventional choice of atomic units h? = 2m = 1 for constant mass m. 
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mass function as well as of the coefficient functions in the ladder operators. We hope to address some 
of these issues elsewhere. 
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